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Abstract. We present here the results of astrometric, photometric and spectroscopic observations leading to the 
determination of the orbit and dynamical masses of the binary L dwarf 2MASSW J0746425+2000321. High 
angular resolution observations spread over almost 4 years and obtained with the Hubble Space Telescope (HST), 
the ESO Very Large Telescope (VLT), and a the W. M. Keck Observatory (Keck) allow us to cover ^36% of the 
period, corresponding to 60% of the orbit, and, for the first time, to derive a precise estimate of the total and 
individual masses of such a late-type object. We find an orbital period of 385O.9I7Q7 days. The corresponding total 
mass is 0.1461q qJ5 Mq, with uncertainties depending on the distance. Spatially resolved low resolution optical 
(550-1025 nm) spectra have been obtained with HST/STIS, allowing us to measure the spectral types of the two 
components (L0=b0.5 for the primary and L1.5=b0.5 for the secondary). We also present precise photometry of the 
individual components measured on the high angular resolution images obtained with HST/ ACS and WFPC2 
(visible), VLT/NACO (J, H and bands) and Keck I (K^ band). These spectral and photometric measurements 
enable us to estimate their effective temperatures and mass ratio, and to place the object accurately in a H-R 
diagram. The binary system is most likely formed by a primary with a mass of 0.085±0.010 Mq and a secondary 
with a mass of 0.066di0.006 Mq, thus clearly substellar, for an age of approximately 300ibl50 Myr. Ha variability 
indicates chromospheric and/or magnetic activity. 

Key words. - stars: very low-mass, brown dwarfs - star: 2MASSW J0746425+2000321 - Binary: orbit, dynamical 
mass 
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gram 70.D-0773), and the W. M. Keck Observatory. 



1. Introduction 

With spectral properties intermediate between those of 
giant planets and late-type stars, ultra-cool and brown 
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dwarfs have opened a new chapter in the study of atmo- 
spheric physics and chemistry. One of the ultimate goals 
of a theory of very low mass and sub-stellar objects is an 
accurate determination of the mass of an object based on 
spectroscopic characteristics and luminosity. The degen- 
eracy in the mass-luminosity (age-temperature) relation 
for ultra-cool dwarfs makes it difficult to pin down their 
physical properties. Luminosities and effective tempera- 
tures of ultra-cool dwarfs are function of both age and 
mass ([Burrows et al.l[l997l) so that an older, slightly more 
massive ultra-cool dwarf can exhibit the same effective 
temperature as a younger, less massive one. Therefore, dy- 
namical masses, which are model-independent, are highly 
required in order to calibrate the mass-luminosity rela- 
tion. The study of binary ultra-cool dwarfs offers several 
advantages for such a study. Both components of the bi- 
nary are expected to be coeval, thus removing part of the 
above mentioned degeneracy in the mass-luminosity (age- 
temperature) relation. Only very few observational con- 
strains on the masses of this class of objects are avail- 
able nowadays, and we present here the first measure- 
ment for field L-dwarfs, at the stellar/sub-stellar tran- 
sition. The only one similar measurement available up 
to now concerned the M8.5/ M9.Q brown dwarf binary 
GJ 569Bab, presented bv Lane et"al] ()200lh . Although the 
age of 2MASSW J0746425-h2000321 is not known inden- 
pendently from any models yet, these observations give 
promising results as a first step toward the calibration of 
the models. 

In this paper, we will present the results of high 
angular resolution imaging and low resolution spectro- 
scopic observations of a binary L dwarf, leading to 
the determination of its orbital parameters, and total 
and individual masses. In section [SJ we will present 
2MASSW J0746425+2000321, in sectionOlwe wih describe 
the observations and the processing of the data, in section 
^we will explain the calculation of the orbital parameters 
and in section [^1 we will present the analysis of the indi- 
vidual spectra and luminosities and discuss the results. 

2. 2MASSW J0746425+2000321 

2MASSW J0746425+2000321 has been observed and re- 
ported in several catalogues a nd articles. It h a s bee n 

1o^, 



identified as a L0.5 dwarf by IKirkpatrick et al 

ll 



and suggested to be a binary by i Reid et al 




based on its position in a colour-magnitude diagram. It 
has b een resolved as a multiple system by Reid et al. 

with a separation of 0'.'22 and a position an- 
gle (P. A) of 15° , a measurement later corrected by 
iBouv et a1.' (2003) to 0'.'219±0'.'003 and P.A=168?8±0?3. 
2MASSW J0746425+2000321 has b een reported insev- 
eral surveys, such as USNO-B (iMonet et al.1 l2003l) . 
GSC2.2, and 2MASS (Cutri et al."2003). Table □ gives an 
overview of the astrometric and photometric properties of 
2MAS SW J0746425 -h2000321 as reported in these cata- 
logues. iDahn et al.1 (^2002) and the USNO-B. 1 catalogue 
both report a measurement of the proper motion of this 



objects, with /iQ,=-370±4 mas yr and /i(5=-42±4 mas 
yr-i fUSNO-B.l) and /i^=-374zbQ 3 mas yr-^ and fis=- 
58±0.3 mas yr" -1 (jPahn et al."2002V These measurements 
confirm that 2MASSW J0746425+2000321 is a common 
proper motion pair. Such a proper motion indeed implies 
a motion of ^1% during the 4 years we made the fol- 
low up observations, whereas the separation between the 
two components varied only of ^O^.^l. Using high reso- 
lution spectra obtained at Keck. IPeid et al.l |2002l ) mea- 
sured a rotational veloci ty of 24 km/s . Using VLT/UVES 
high resolution spectra. iBailer- JonesI (accepted for A&A, 
I2OO4I) measured a rotational velocity ranging between 
25. 6< vsini <30.6 km/s, corresp onding to a period 
bet ween 1.73<T<3.71 h ours. Both IClarke et all ()2002l) 
and Gelino et al. (2002) report photometric variability, 
which they attribute mainly to the for mation of cl o uds in 
the upper layers of the atmospheres. IPahn et al.l ()2002l ) 
measured its distance using trigonometric parallax at 
12.21±0.05 pc. 



3. Observation and data processing 

Table El gives a log of all the observations we use in this 
study. 



3.1. High Angular Resolution imaging with H ST/ ACS 
and STIS 

High angular resolution images have been obtained 
with the HST Advanced Camera for Survey (ACS, 
IPavlovskv et al.l UoO."^) and Space Tele scope Imaging 
Spectrograph fSTIS. Ixim Ouiiano et "Zl Eo03). We ob- 
served 2MASSW J0746425+2000321 using the ACS and 
its High Resolution Channel (HRC) in three different opti- 
cal filters (F625W, F775W and F850LP), and STIS in the 
long-pass filter. 2MASSW J0746425+2000321 is clearly re- 
solved on both sets of data (see Figure^]), and we were able 
to get precise astrometric and photometric measurements. 
The data have been analyzed using a custom-made pro- 
gram performing PSF fitting to compute the precise sep- 
aration, position angle and flux ratios of the multiple sys- 
tem. T his program is identical to that used bv'B ouv et al.l 
()2003l) . adapted for ACS/HRC and STIS, and is fully de- 
scribed in the mentioned paper. A brief summary is given 
here for completeness: the PSF fitting routine builds a 
model binary using ten different PSF stars coming from 
different ACS/HRC images. A cross-correlation between 
the model and the binary system yields the best values for 
the free parameters: separation, position angle and fluxes. 
The uncertainties a nd limitations of th is technique are dis- 
cussed in detail in 'Bouv et al.l (l2003[). The pr ogram has 
been slightly improved since B ouv et al.l ()2003l ) presented 
it in their paper, and the systematic errors and uncertain- 
ties on the difference of magnitude are now much better 
understood and corrected (Bouv 2004). While the original 
program was using only 5 free parameters (the coordinates 
of the primary, the coordinates of the secondary, and the 
flux ratio between the two components), the new version 
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uses 6 free parameters (the 4 coordinates, the flux ratio 
between the primary and the PSF star, the flux ratio be- 
tween the secondary and the PSF star), which allows to 
decrease considerably the systematic errors and uncertain- 
ties. 

3.2. High Angular Resolution imaging with 
VLT/NACO 

We also obtained high angular resolution images using 
the ground based facilities offered by ESO on Cerro 
Paranal on 2003 February 18th and 2003 March 22nd. The 
VLT on Yepun uses NACO, an adaptive optics platform 
(IPousset et a1.lE003[ iLenzen et a].ll2003t iBrandner et a1.l 
l2002j) to achieve diffraction limited images. NACO offers 
the possibility to use an infrared wavefront sensor, and 
is therefore ideally suited for the study of ultra-cool and 
red objects. Its CONICA array offers a 0'.'01326±0'.'001 
pixel-scale that provides critical Nyquist sampling of the 
diffraction limited images of the telescope at these wave- 
lengths. Its absolute orientation is known to within 

The atmospheric conditions during the observations 
were good (respectively A/ro=0'.'62 and airmass=1.5, 
and X/ro=0%7 and airmass=1.4), and very sharp im- 
ages in (first observation) and J, H and (sec- 
ond observation) were obtained with strehl ratios of 
Sr{Ks) -30% (first obs.) and Sr{J) -13% ,Sr{H) -27%, 
and Sr{Ks) —46% (last obs.). Figure H shows the two 
images obtained during these two nights. 

During the last observation, a PSF star was also ac- 
quired in order to perform accurate photometry of the 
adaptive optics data of the corresponding night. The 
object, DENIS-P J131500.9-251302 (spectral type -M8, 
J=15.2, 11=14.54 and K5=14.02 mag), was observed 
under better conditions (A/ro=0'.'43 and airmass=1.03) 
with a strehl ratio of Sr{H) -10%, and Sr{Ks) -40%. 
Unfortunately it was not observed in J. We performed 
the photometry using standard DAOPHOT PSF fitting 
photometry. The results are summarized in Table [2l 
Although the PSF star has a spectral type earlier than 
2MASSW J0746425+2000321 and was observed at much 
better airmass and better seeing, the relative photom- 
etry we obtain in H and is in very good ag:ree- 
ment with the one reported by IClose et al.l ()2003l) with 
Gemini North/Hokupa'a and the one we measure with 
Keck I/NIRC {Ks band). 

3.3. Speckle Observations with Keck 

On 2003 December 04, we obtained K band speckle obser- 
vations of our target at the 10 m Keck I telescope w ith the 
facility instrument NIRC (iKleinmann et al]ll994^ ■ With 
its re-imaging optics ('Matthews et al."l996^, this 256x256 
near-infrared array offers a 0'.'0203±0'.'0003 pixel scale that 
provides Nyquist sampling of the diffraction limit of the 
telescope at this wavelength (about 0'.'05); its absolute 
orientation is known to within 1°. Several stacks of 200 



short-integration exposures were obtained (t — 0.1 s, i.e., 
fast enough to effectively "freeze" the atmospheric turbu- 
lence and retain the high-angular resolution information 
in the image), and similar stacks on two calibration point 
sources were obtained immediately before and after our 
target. Standard speckle data reduction routines wer e ap- 
plied to the data; we refer the reader to iGhej ()l993l) and 
IPatience et al.l ()l998l) for more details and only summarize 
briefly the various stages involved in the data reduction 
process. Each individual exposure is first sky subtracted, 
flat-fielded and bad pixel-corrected; its power spectrum is 
then calculated. The power spectra are median- averaged 
over each stack and divided by that of the calibrator. A 2- 
D sinusoidal function is then fitted to the power spectrum 
to determine the binary properties: separation, position 
angle and flux ratios. Uncertainties are estimated from 
the standard deviation of the parameters extracted from 
all stacks. There is a 180° ambiguity in the position angle 
of the binary as derived through power spectrum analysis, 
but this can be resolved by shift-and-adding all individ- 
ual exposures using the brightest speckle as a reference. 
The resulting image shows the companion to be roughly 
to the Northeast of the primary and the astrometric ac- 
curacy of the power spectrum analysis is much higher. 
The results obtained are reported in Tables 01 and El and 
Figure m shows the final image. As mentioned above, the 
difference of magnitude is in very good agreement with 
the previous measurements within the uncertainties. 

3.4. High Angular Resolution - Low Spectral 
Resolution Optical Spectroscopy with STIS 

In order to get the spatially resolved spectra of each com- 
ponent of the multiple system, we tried to align the slit 
along the axis of the binary. Scheduling constraints of 
HST made it difficult to get long slit STIS observations of 
2MASSW J0746425+2000321 at a particular roh angle of 
HST. In order to ease scheduling, a range of admitable roll 
angle was defined. This combined with the relatively rapid 
orbital motion of 2MASSW J0746425+2000321 meant 
that the entrance slit of STIS was not optimally aligned 
with the position angle of 2MASSW J0746425+2000321 
at the time of the observations. Since the size of the slit 
we used (0'.'2) is larger than the separation of the binary 
(— 0'.'125), we could nevertheless obtain a resolved 2-D 
spectrum and perform an extraction of the two spectra. 
The effects on the spectral analysis can be the following: 
since the red and the blue photocentres of each component 
are not symmetrically centred in the slit, the dispersion 
of the red and blue parts of the spectrum suffers differ- 
ently from obstruction by the slit. Since the position of 
the photocentres and the dispersion of the light depends 
on the wavelength, the loss in flux also varies with the 
wavelength. This effect produces a "bluer" spectrum for 
the secondary. 

The separation between the two spectra is about 2 pix- 
els, whereas the full width at half maximum (FWHM) of 
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the line spread function (LSF) varies between 1.0 and 1.2 
pixel, so that the two spectra are barely resolved. In order 
to extract the spectrum of each component, we used a cus- 
tom made program able to perform a fit of the two blended 
spectra. On each cross-dispersion column, a minimum 
fit was performed to the data using the cross-dispersion 
profile of a reference spectrum at the same wavelength. 
The latter spectrum was obtained with the same instru- 
ment settings on a K7 dwarf (TWA 6, program 8176, P.I. 
Schneider). The free parameters for the fit are the ampli- 
tude of the primary, the amplitude of the secondary, the 
position of the primary and the position of the secondary. 
Since the cross dispersion profile is barely sampled, we 
also performed a linear re-sampling of the data by a fac- 
tor of eight prior to the fit, in order to avoid problems 
due to spectral aliasing. To ensure more robustness and 
increase reliability, the program was used in 2 passes. A 
polynomial fit of the results on the positions of the two in- 
dividual spectra was made after the first pass, in order to 
identify and remove outliers (due to bad pixels or cosmic 
rays). The results of these fits were then used as first guess 
inputs for the second pass. The results obtained with the 
second pass are very close to that obtained with the first 
pass but cleaner (without the more obvious bad pixels, 
cosmic rays and outliers), ensuring that the whole algo- 
rithm is robust enough and converging properly. Figure [2l 
gives an overview of the results at a particular wavelength. 
The residuals after the fit represent between 1.5% and 9% 
of the total intensity of the original spectrum depending 
on the wavelength, thus of the same order than the signal 
to noise ratio, ensuring that the quality of the fit is good. 

4. Orbital parameters and determination of the 
total mass 

We used three different and independent custom-made 
programs to determine the best fitting orbital solution 
for 2MASSW J0746425+2000321 and the uncertainties on 
each of the fitted parameters. The orbit can be entirely 
described by seven independent parameters: semi- major 
axis (a), orbital period (P), inclination (i), eccentricity 
(e), position angle of the ascending node (f^), angle be- 
tween the ascending node and periastron {uj) and time 
of periastron passage (Tq). With seven two-dimension as- 
trometric data-points, this fully-constrained problem has 
seven free parameters. The total mass of the system can 
be derived from the orbital period and semi-major axis 
through Kepler's Third Law. 

4.1. ''Amoeba" method 

The first method minimizes in the nonlinear 7-dimensions 
function by down hill simplex metho d, using the amoeba 
algorithm (see e.g 'Pres s et al.l[l992l for a description of 
the method and algorithm.). It fits all seven orbital param- 
eters simultaneously, taking into account non-equal errors 
of the measurements. The reduced-x^ of 1.41 ensures that 
the fit is satisfactory (see section lOI below) . The results 



are shown in Figure and Table [SI No uncertainties on 
the derived parameters are available with this method. 

4.2. Iterative method 

This method uses 50 000 independent starting points that 
consist of a set of the 7 parameters being randomly cho- 
sen from their entire range of possible values. For each 
starting point, a Powell converge nce algorithm minimizing 
the total (|Press et al.lll993) modifies simultaneously 
all 7 parameters until it converges to a local minimum. 
Once convergence for all 50 000 sets of initial guesses has 
been achieved, we read through the output file to find the 
absolute minimum of the function, which reveals the 
best-fitting orbital solution. Our best-fitting solution, il- 
lustrated in Figure ini has a satisfying reduced-x^ value of 
1.38. 

Uncertainties for each parameters are defined by the 
range of possible values indicated by all solutions with to- 
tal x^ between Xmin Xmin + 4- These represent the 
95.4% confidence level for each parameters. Due to the 
highly non-linear behavior of the equations of orbital mo- 
tion, the uncertainties do not follow a Gaussian statis- 
tics and nor even symmetric about the best fit. Note 
that the uncertainties are derived under the assumption 
that all parameters are independent, which is not cor- 
rect. For instance, the uncertainties derived for P and a 
would yield an uncertainty on the system mass on order 
of ±0.055 M0, ~4 to 9 times larger than we actually de- 
rived here. Therefore, the uncertainties quoted here are 
only valid if they are used for one parameter at a time. 
Figure ^ shows that this is because the fitted values of 
P and a are tightly correlated and correspond to a very 
narrow range of possible masses. Although the orbit is 
not perfectly known yet, the total mass is relatively pre- 
cisely determined: Mtot=0.146to.oo6 corresponding 
to 4~11% uncertainty. 

The uncertainty on the distance to the system trans- 
lates into a separate 2.4% uncertainty (2-cr), or 0.035 Mq, 
on the system mass. The uncertainty on the orbital fitting 
is therefore the major source of uncertainty for this binary 
system since both sources of uncertainty (fit and distance) 
should be added in quadrature. 

4.3. ORBIT 

We a lso used the ORBIT program of iForveille et al.l 
(1999), fully described in their article. Briefly, "t/ie pro- 
gram performs a least square adjustment to all available 
observations, with weights inversely proportional to the 
square of their standard errors [. . .] The program uses a 
Levenberg-Marquardt minimization algorithm i Marauardl 
\l96i) [...] Standard errors for derived parameters are com- 
puted from the full covariance matrix of the linearized least 
square adjustment.''^ 
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4.4. Reduced-x^ and uncertainties 

The reduced-x^ values of ~1.4 indicate that some of the 
uncertainties on the astrometric measurements may be 
slightly underestimated, although with only seven mea- 
surements for seven free parameters, such a value is sta- 
tistically acceptable. Although rescaling the astrometric 
uncertainties to reach a reduced-x^ of 1.0 could be argued 
for, the diversity of the instruments used for this orbital 
analysis suggests that such a treatment would be at least 
as erroneous as it could be helpful. For the time being, we 
decided to stick to the quoted astrometric uncertainties to 
derive the uncertainties on the orbital parameters of the 
binary. 

5. Discussion 

5.1. Spectral Types, Effective Temperatures 

The composite spectrum of 2MASSW J0746425- 
+2 000321 AB in the optica l has been previously studied 
bv IKirkpatrick et al.l ()2000l) . who derived a spectral type 
of L0.5. 

In order to derive the spectral types of each compo- 
nent, we compared their STIS spectra, extracted with 
the procedure described in Section 13.41 with the spec - 
tra of field L dwarfs published by Ma rtin et al.l ()l999l) . 
As explained in this latter article, the relative strength 
of the TiO bands between 840-860 nm and the CrH 
and FeH bands between 860-880 nm are good indica- 
tors for the effective temperature changes. In later L 
dwarfs the TiO bands get weaker with respect to the 
CrH and FeH bands. Figures [SI and El shows the com- 
parison. 2MASSW J0746425+2000321 A is clearly very 
similar to the LO field dwarf DENIS-P J090957.1-065806, 
while 2MASSW J0746425+2000321 B is between the 
LI (DENIS-P J144137.3-094559) and L2 (Kelu 1) field 
dwarfs. We thus derive a spectral type of LO+0.5 and 
LI. 5 ±0.5 for A and B, respectiv ely. This is consistent wit h 
the spectral type obtained bv IKirk patrick et al.' (2000), 
which is a blend of A and B. It is also consistent with the 
modest difference in brightness (see Table |2|, which im- 
plies a difference in temper ature of o nly 100 K according 
to the models of,Chabrier et"Z| ^0(\). 

5.2. Spectral Features 

The main emission lines present in the primary's spectrum 
is Ha (EW=-25.0±0.5 A, l-a). The spectrum of the sec- 
ondary is more noisy but the Ha emission line appears 
clearly, with an equivalent width of EW=-18 . 0+0.5 A. 
From their high-resolution spectra, Reid et al.l (|2002) re- 
ported a Ha emission of -1.2 A for the unresolved system. 
The difference between the two measurements indicates 
that 2MASSW J0746425+2000321 A and B display some 
chromospheric and/or magnetic activity. Li I absorption 
is not detected in any of the two comp onents with an 
upper limit of ~1.5 A. lR,eid et all |2002) did not detect 
any Li I absorption with an upper limit of detection at 



<0.5 A. The presence of strong resonance doublets of al- 
kali elements (K I at 766.5 and 769,9 nm; Na I at 818.4 and 
819.5 nm; and Cs I at 852.1 and 894.3 nm) as well as strong 
metallic molecular band- heads (CrH at 861.1 n m and FeH 
at 86 9.2 nm) is characteristic for L-dwarfs (Ma rtm et al.l 
Il997l) . The measurements of equivalent widths of the main 
atomic lines are reported in Table El It is interesting to 
note that the equivalent width we measure for Cs I at 
852.1 nm for the primary corresponds to an effective tem- 
perature o£^^4900-2000 K in the effective temperature 
scale of iBasri et al.l (|2000), therefore in good agreement 
with t he effective temperature derived bv lSchweitzer et all 
(|200ll) from their comparison of low and high resolution 
Keck spectra with the DUSTY models. 

5.3. Colour-Magnitude Diagrams 

Figure [3 shows a colour-magnitude diagram of 2MASSW- 
J0746425+2000321 AB, 2MASSW J0746425+2000321 A 
and 2MASSW J0746425+2000321 B, and compares with 
the isochrones of the most recent DUSTY models for solar 
metallicity. To convert the observed magnitudes to abso- 
lute magn itudes we used th e trigonometric parallax re- 
ported bylDahn et al.| §^1). 

The position of 2MASSW J0746425+2000321A shows 
that the age ranges between 150 and 500 Myr, thus 
relatively young. This is not c onsistent with the ver y 
high surface gravity obtained bv ISchweitzer et "all ()200lL 
They used high and low resolution unresolved spectra 
and compared it w ith the DUSTY atmospheric models of 
lAllard et al.l ()200ll ). Using a x^-fitting algorithm, they ob- 
tain an effective temperature of 2000 K and a surface grav- 
ity log ^=6.0 from their low resolution spectra. From their 
high resolution spectra, they obtain an effective tempera- 
ture of 1900^2000 K and a surface gravity log ^=5.0^5.5. 
The temperatures are in good agreement with the spectral 
types and colours we report here, but the surface gravity 
is too high for the young age. This is probably because 
the DUSTY models overestimate the dust effects. As a 
consequence, the strength of the alkali lines in the optical 
decreases, and the surface gravity is biased toward higher 
values. 

The mass of the primary ranges between 0.075 and 
0.095 M0, while the mass of the secondary has large error 
bars and ranges between 0.055 and 0.100 Mq. The total 
mass of the system therefore ranges between ~0.130 and 
~0.190 Mq, which is consistent with the dynamical mass 
considering the large uncertainties in the H-R diagram. 

For the secondary's mass it appears more appropri- 
ate to use the mass of the primary from the H-R dia- 
gram, which has reasonable uncertainties, together with 
the very precise dynamical total mass. This yields a mass 
between 0.052 and 0.072 Mq, therefore clearly substellar. 
The absence of lithium absorption in the spectra gives also 
a constrain on the lower limit of the mass. According to 
the DUSTY evolutionary models. Lithium should be de- 
pleted for masses greater than 0.075 Mq at 150 Myr and 
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raasses greater than 0.060 Mq at 500 Myr. The mass of 
the secondary must therefore be greater than 0.060, and 
ranges between 0.060<Mi5 <0.072 Mq. 

The system is thus very likely made of a brown dwarf 
orbiting a slightly more massive very low mass star. Both 
objects are very close to the stellar-substellar boundary. 

6. Conclusion 

We present new astrometric, photometric and spectro- 
scopic observations of the ultra-cool 2MASSW J0746425- 
+200032 lAB pair, which have enabled us to compute 
the orbital parameters, total and individual masses and 
spectral types of the system. We find a total mass of 
0- 146^0.006 ^0 with a L0±0.5 primary and a 

L1.5±0.5 secondary. The orbit is eccentric (e = 0.411q q9) 
with a period of 3850.9^767 days. Our observations en- 
able us to follow 60% of the orbit, corresponding to 36% 
of the period. Near-infrared photometry of the individ- 
ual components enables us to locate them in a H-R dia- 
gram and compare them to the most recent evolutionary 
DUSTY models. The pair is likely formed by a brown 
dwarf (L1.5±0.5, with 0.060<Mi5 <0.072 Mq) orbit- 
ing a L0±0.5 very low mass star (0.075<M<0.095 Mq). 
The system appears to be young, with an age in the 
range 150^500 Myr. Further high-angular resolution im- 
ages should allow to cover 100% of the orbit and refine 
the orbital parameters. Radial velocity measurements, if 
feasible, would then allow to compute precise individual 
masses. 
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Table 1. Astrometry and photometry from different catalogues 



Date 


R.A 


Dec. 


Uncert. 


Filter 


Mag. 


Source 


Ident. 


dd/mm/yyyy 


(J2000) 


(J2000) 






[mag] 






01/01/1984^ 


07 46 42.5 


+20 00 32.6 


±ori 


Rl 


18.28 


USNO-Bl.O 


USNO-Bl.O 1100-0150847 










B2 


21.7 














R2 


17.87 






01/01/1998 


07 46 42.55 


+20 00 32.14 


±0^3 


R 


17.6 


GSC2.21 


GSC 2W 22110125398 


05/12/1997 


07 46 42.56 


+20 00 32.2 


±ori 


J 


11.759 


2MASS 


2MASSW J07464256+2000321 










H 


11.007 














K 


10.468 






06/12/2002 








L' 


11.19 


Leggett et al. (2QQ2) 





^ Mean epoch of observation 



Table 2. Observation log. 



Imaging 



Instrument 


Filter 


Exp. 


Date Obs. 


Pixel Scale 






Time [s] 


dd/mm/yyyy 


n 


HST/WFPC2-PC 


F814W 


50 


25/04/2000 


0:'0455 


Gemini North/Hokupa 


J 


120 


07/02/2002 


0'.'0199 


Gemini North/Hokupa 


H 


720 


07/02/2002 


0^0199 


Gemini North/Hokupa 


K 


120 


07/02/2002 


0^0199 


HST/ACS-HRC 


F625W 


960 


21/10/2002 


0^0250^ 


HST/ACS-HRC 


F775W 


440 


21/10/2002 


0^.^0250^ 


HST/ACS-HRC 


F850LP 


340 


21/10/2002 


0^0250^ 


VLT/NACO 


Ks 


0.4 


18/02/2003 


0'.'0148^ 


VLT/NACO 


J 


10 


22/03/2003 


0'.'0148^ 


VLT/NACO 


H 


5 


22/03/2003 


0'.'0148^ 


VLT/NACO 


Ks 


5 


22/03/2003 


0'.'0148^ 


Keck I/NIRC 


Ks 


20 


04/12/2003 


0^.^0203 


HST/STIS 


longpass 


10 


09/01/2004 


0^0508 


Spectroscopy 


Instrument 


Wavelength 


Exp. 


Date Obs. 


Dispersion 




Range [nm] 


Time [s] 


dd/mm/yyyy 


[A/pixel] 


HST/STIS 


525-1300 


1980 


09/01/2004 


4.92 



^ Effective value on the processed images, slightly different from the 0^.^028x0^025 given in the manual. 
^ Effective value on the processed images, slightly different from the 0^.^01326 given in the manual. 



Press, W. H., Teukolsky, S. A., Vetterling, W. T., & 
Flannery, B. P. 1992, Numerical recipes in C. The art 
of scientific computing (Cambridge: University Press, 
— cl992, 2nd ed.) 

Reid, I. N., Gizis, J. E., Kirkpatrick, J. D., & Koerner, 
D. W. 2001, AJ, 121, 489 
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AJ, 119, 369 

Reid, 1. N., Kirkpatrick, J. D., Liebert, J., et al. 2002, AJ, 
124, 519 
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in Adaptive Optical System Technologies 11. Edited 
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Table 3. Relative Photometry of 2MASSW J0746425+2000321AB 



Date 


Instrument 


Filter 


Mag. Prim. 


AMag. 


Source 


dd/mm/yyyy 








[mag] 


[mag] 




25/04/2000 


HST/WFPC2 


F814W 


15.41=b0.15 


1.00±0.09 


(1) 


07/02/2002 


Gemini North/Hokupa'a 


J 


12.19±0.07 


0.60±0.20 


(2) 


07/02/2002 


Gemini North/Hokupa'a 


H 


11.54=b0.11 


0.48=i=0.15 


(2) 


07/02/2002 


Gemini North/Hokupa'a 


K' 


11.05±0.09 


0.44±0.15 


(2) 


21/10/2002 


HST/ACS 


F625W 


18.81=t0.05 


0.48±0.03 


(3) 


21/10/2002 


HST/ACS 


F775W 


15.98±0.05 


0.68±0.04 


(3) 


21/10/2002 


HST/ACS 


F850LP 


14.24±0.05 


0.76±0.04 


(3) 


22/03/2003 


VLT/NACO 


H 


11.55±0.08 


0.46±0.15 


(3) 


22/03/2003 


VLT/NACO 


Ks 


11.06±0.09 


0.42±0.15 


(3) 


04/12/2003 


Keck I/NIRC 


Ks 


11.03=b0.03 


0.52=b0.03 


(3) 


^ Source: m IBouv et alJ (1200^1: r2)IClose et all tOO^: (S) this Dauer 






Table 4. Relative Astrometry of 2MASSW J0746425+2000321AB 




Date 


Time^ 


Sep.2 


PA^ 


Instrument 


Source^ 


dd/mm/yyyy 


hh/mm/ss 


[mas] 


n 








25/04/2000 


08:14:14 


219±3 


168.8=t0.8 


HST/WFPC2 


(1) 


07/02/2002 


09:48:55 


121zb8 


85.7±3.5 


Gemini North/Hokupa'a 


(2) 


21/10/2002 


23:10:43 


119.5±1 


33.9±0.5 


HST/ACS 


(3) 


18/02/2003 


01:40:45 


131.3=1=3.9 


13.8=bl.9 


VLT/NACO 


(3) 


22/03/2003 


01:22:00 


123.5=b2.1 


4.6=bl.O 


VLT/NACO 


(3) 


04/12/2003 


15:15 


126.5=bl.8 


317.9=b0.7 


Keck I/NIRC 


(3) 


09/01/2004 


18:51:45 


134.5±3 


311.1±1.2 


HST/STIS 


(3) 



^ The uncertainty corresponds to the exposure time (see Table |2j 

^ 1-cr uncert ainties (combined i nstru mental and measu rement) 

^ Source: m IBouv et Zl (|20n,1) : r2) IClose et Zl(l200,1) : (3) this paper 



by Wizinowich, Peter L.; Bonaccini, Domenico. 
Proceedings of the SPIE, Volume 4839, pp. 140-149 
(2003)., 140-149 
Schweitzer, A., Gizis, J. E., Hauschildt, P. H., Allard, F., 
& Reid, 1. N. 2001, ApJ, 555, 368 
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Table 5. Orbital Parameters of 2MASSW J0746425+2000321AB 



9 



Parameter Iterative Method^ Amoeba Method ORBIT'^ 



Total Mass [M^t)] 




0.1511 


0.148 


Period, P [days] 


3850.91?°^ 


3718 


3863±609 


Eccentricity, e 


o.4il°:°9 


0.3999 


0.417±0.062 


Semi-major axis, a, [A.U] 


9 Kq+0.37 
■^•^^ — 0.28 


2.50 


2.55±0.25 


Inclination, i [°] 




141.65 


140.65±2.29 


Argument of Periapsis, u [°] 




350.15 


350.65±3.58 


Longitude of ascending node, Q [°] 




18.97 


20.84=b7.68 


Periastron Passage, To (year) 


2002.89l°;J^ 


2002.91 


2002.84zb0.07 


reduced-x^ of the fit 


1.38 


1.41 


1.46 



^ 2-cr uncertainties, corresponding to a 95.4% level of confidence. These uncertainties do not include the uncertainty on the 
distance (~2.4%, 2-cr), which should be added in quadrature. 

^ 2-cr uncertainties. These uncertainties do include the uncertainty on the distance, but assume that the uncertainty are linear, 
which is not the case here. 



Table 6. Atomic lines in the spectra of 2MASSW J0746425+2000321A and B 





K I Wavelength 


Cs I Wavelength 


Na I D Wavelength^ 


Object 


7665 7699 


8521 8943 


8183-8195 


2MASSW J0746425+2000321A 


22.6 17.4 


2.1 0.99 


8 


2MASSW J0746425+2000321B 


19.0 16.4 




5 



Note. — All units are in angstroms. 1-cr uncertainties are ^^0.5 A. 
^ Corresponds to the blend of the 8183 and 8195 A doublet. 
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Fig. 1. Images of 2MASSW J0746425+2000321A and B obtained at different epochs with HST, Gemini, NACO and 
Keck 1. The scale and the orientation are the same for all images, and indicated on the figure. 
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Fig. 2. Extraction of the individual spectra. This figure shows the cross dispersion profile of the 2-D spectrum around 
850 nm, and the best fits of the primary and the secondary (light blue). The sum of the primary and secondary (red) 
is also indicated for comparison with the raw data (black). In that case the intensity of the residuals (green) is less 
than ~2% of the intensity of the raw data. 
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Fig. 3. Positions of 2MASSW J0746425+2000321A and B and best fit of the orbit . The dotted curve represents the 
best bit orbit obtained with the amoeba method, and the sohd curve the result obtained with the iterative method, 
and the dashed curve the solution given by ORBIT. It appears clearly that the three methods give close results, well 
within the uncertainties. The plus indicate the observations and their uncertainties, and the corresponding epoch is 
indicated. The central cross shows the position of the primary. 
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Orbital Period (yrs) 

Fig. 4. Ambiguity in tlie orbital parameters: tlie contours represent tlie solutions at 1- and 2-cr (Ax^=1.0 and 4.0 

respectively) in the semi-major axis vs period space. The lines corresponding to different masses indicate the levels 
of confidence at 68.7% and 95.4% (1- and 2-(7 respectively). The filled square in the centre corresponds to the best 
fit. Although relatively large ranges are possible for the period and the semi-major axis, the range of corresponding 
acceptable masses is very narrow: between 0.143 and 0.153 Mq at 1-a (represented by dashed lines) and between 0.141 
and 0.160 Mq at 2-a (dotted lines). Although the orbit is not perfectly known yet, the mass is relatively precisely 
determined, with a best value at 0.146lo.oo6 ^0 (^"^ uncertainties, solid line). 
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Fig. 5. STIS optical low resolution spectra of 2MASSW J0746425+2000321A compared to spectra of filed ultra-cool 
dwarfs. The four plots show the STIS spectrum of 2MASSW J0746425+2000321A, smoothed via a boxcar (width = 
5 pixels), and compared to: a) DENIS-P J024351.0-543219 (dM9); b) DENIS-P J090957.1-065806 (dLO); c) DENIS- 
P J1 44137. 3-094559 (dLl ); d) Kelu 1 (dL2). All spectra have been normalized at 840 nm. Spectra for the field dwarfs 
fromlMartm etZj] 
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Fig. 6. Same as Figure but for 2MASSW J0746425+2000321B 




Fig. 7. Colour-Magnitude diagrams M^^ vs (J - Ks) displaying the location of 2MASSW J0746425+2000321A, B 



